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ABSTRACT 

The asymmetric molecular emission lines from dense cores reveal slow, inward mo- 
tion in the clouds' outer regions. This motion is present both before and after the 
formation of a central star. Motivated by these observations, we revisit the classic 
problem of steady, spherical accretion of gas onto a gravitating point mass, but now 
include self-gravity of the gas and impose a finite, subsonic velocity as the outer bound- 
ary condition. We find that the accretion rate onto the protostar is lower than values 
obtained for isolated, collapsing clouds, by a factor that is the Mach number of the 
outer flow. Moreover, the region of infall surrounding the protostar spreads out more 
slowly, at a speed close to the subsonic, incoming velocity. Our calculation, while highly 
idealized, provides insight into two longstanding problems - the surprisingly low accre- 
tion luminosities of even the most deeply embedded stellar sources, and the failure so 
far to detect spatially extended, supersonic infall within their parent dense cores. In- 
deed, the observed subsonic contraction in the outer regions of dense cores following 
star formation appears to rule out a purely hydrodynamic origin for these clouds. 

Subject headings: ISM: clouds, kinematics and dynamics — stars: formation — accretion 



1. Introduction 

Astronomers are elucidating, in much greater detail than ever before, the earliest phases of 
low-mass star formation. It has long b een known that at leas t half of dense cores within molecular 



clouds do not yet contain young stars (IBeichman et al.l ll986). It is also known that these starles s 



cores have, on average, less mass than cores with embedded stars (jJijina. Myers. &: Adamsl ll999). 
What is now clearer is the chemical and structural progression of cores as they approach the col- 
lapse threshold and then cross it. The most centrally concentrated starless cores tend to exhibit 
the asymmetric emission line profile s signifying inward contracti on, with the inferred velocity being 



several tenths of the sound speed (|Lee. Myers, fc Tafallal l200ll ). The same blueward asymmetry 
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is seen in those cores with embedded stars, where it is especially common (jMardones et al.lll997l ) 
Detailed examination of line profiles reveals no qualitative change across the collap se threshold, ex- 



cept for the appearance of broad wings that are attributable to molecular outflows (jGregersen et al 
20001 ) . 



It is likely that the contraction observed in starless cores is driven by sel f-gravity, as material 
in the core's outer region settles toward the center (e.g., ICiolek fc Basull200ll ). A similar, exterior 
contraction should also occur in cores that have produced stars, and observations thus far support 
this conclusion. Here, targeted studies a re rare, in p a rt bec ause outflows often contaminate the 
emission line spectra. As one example, iMvers et al.l (|19961 ) utilized a simple radiative transfer 
model to measure the spatially averaged contraction speed in L1526, a core with an embedded star 
that drives a weak outflow. Their derived contraction speed of 0.025 km s -1 is indeed close to 
the inferred speeds in starless cores . The CS and H2CO spectra for the low-mass infall candidates 
identified by iMardones et al.l (|1997l ) have a similar degree of asymmetry as L1526, but have not 
been modeled to this level of detail. 

Taken as a whole, these observations suggest that a dense core reaches the point of collapse 
by accruing matter externally. The core its elf is a relatively quiescent st ructure, since its observed 
emission lines have nearly thermal widths (IBarranco fc Goodman! Il998l ). The exterior molecular 
gas, however, has typic al speeds of 1 to 2 km s _1 , corresponding to Mach numbers of 5 to 10 
(IBergin fe Tafallal 120071 ) . How exactly material passes from this envir onment into the co re is un- 
known, despite increasingly detailed mapping of the transition region (jPineda et al.ll2010l ). In any 
event, the core is dynamically evolving while it collects mass. Gas in its outer mantle creeps inward, 
and continues to do so even after the central region collapses to form a star. 

The standard theory of collapse, first formulated prior to the disc overy of de nse cores, paints 
a very different picture. In the still widely used self-similar model of IShul (j 19771 ). the cloud is a 
perfectly static entity before the first appearance of a protostar. The impulsive suctioning of gas 
caused by this event sends out a sharp rarefaction wave, expanding at the speed of sound. It is 
within this wavefront that infall occurs. Gas leaves the rarefaction wave with zero velocity and 
attains the sound speed roughly halfway inside it. 

For a typical dense core radius of 0.1 pc and sound speed of 0.2 km s -1 , the rarefaction 
wave reaches the boundary in 0.5 Myr. This interval, or one comparable, has been traditionally 
considered to be the lifetime of the protostar phase. Comparing the dynamical picture from theory 
with observations, it is striking that there are no known cores with embedded stars in which the sonic 
transition occurs reasonably close to the cloud boundary. On the contrary, the inferred location of 
this point in the best spectroscopic collapse candidates is only 0.0 1 -0.02 pc from the pr otostar itself 



( Choi et al. 



1995 



Gregersen et al.lll997l : |Pi Francesco et al. 



2001 



Belloche et al.ll2002l ). Again, we 



must remember that strong molecular outflows can both mask and mimic large-scale motion in 
these studies. On the other hand, submillimeter continuum ma pping of cores with embedded star s 
also indicates that infall is confined to a small, central region. IShirlev. Evans. fe Rawlingd ([2002) 
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have shown that more widespread collapse would give rise to a change in the interior density profile 
that is not observed. The continuing absence of dense cores that contain stars and exhibit true 
global infall is forcing us to conclude either that all observed protostars happen to be especially 
young, or else that infall actually spreads more slowly than traditional theory posits. 



The collapse model of lShul (|1977l ) also predicts that the mass accretion rate onto the protostar 
is a constant in time: 

c 3 

M = m -^ . (1) 



Here, c s is the isothermal sound speed in the core and m a is a nondimensional number with a value 
of 0.975. Many authors, starting with lKenyon et al.1 (|1990l ). have noted that the accretion luminos- 
ity associated with this rate exceeds those observed for embedded, infrared sources. Subsequent 
numerical simulations relaxed the assumption of self-similarity, but usually imp osed a rigid bound- 
ary on the dense core, i.e., a fixed surface where the fluid velocity v anishes (e.g. jFoster h Chevalier 



1993 



Ogino. Tomisaka. Nakamuralll999 ; Vorobvov h Basull2005 ). In these calculations, M varies 



in time, but is even higher, at least before gas begins to drain fro m the boundary region. T he sonic 
transition occurs only slightly inside this boundary (see Fig. 1 of lFoster fc Chevalierlll993l ). 



In response to the persistent and growing disparity between theory and observation, researchers 
have taken two very different paths. Some have focused on the possibility that, while cloud collapse 
occurs in the traditional manner, infalling gas is retained in a circumstellar d isk, and only released 
periodically in bursts onto the protostar (e.g., iDunham &: Vorobvovl 120121 ). If the duration of 
bursts is sufficiently brief, we could be viewing protostars only when M is relatively low. Whatever 
the merits of this idea, it does not address the troubling aspects of collapse models generally. 
Protostellar luminosity and cloud collapse are clearly linked, and a solution that considers both is 
preferable. 



Fo l lowing this second, unified path, iFatuzzo. Adams. &: Mversl (|2004l ) generalized the model 
oflShul (|1977l ) and i nvestigated th e self- simila r evolution of a c l oud t hat has a specified velocity 
profile at the start. iGomez et al.l (|2007l ) and iGong Qstrikerl (|2009l ) simulated numerically the 
birth of a dense core as arising from the co nvergence of a s upers onic inflow. Both groups tracked 
the evolution past the start of core collapse. IStahler Yen! (|2009l ) used perturbation theory to find 
the contraction in starless cores that are just ap proaching collapse, and were able to reproduce the 
asymmetric emission lines (jStahler &: Yenll2010l ). 



In this brief contribution, we focus on how continuous, external mass addition to the dense 
core affects protostellar infall itself. We show that M in this case is equal to c 3 /G times a factor 
that is essentially the Mach number of the incoming flow. This key result is easy to understand 
physically. For a cloud that is marginally stable against self-gravity, the free-fall and sound-crossing 
times are comparable. It follows that the mean speed of infalling fluid elements in an isolated cloud 
undergoing collapse is approximately c s . When accretion is imposed from the outside, this mean 
speed is closer to the externally impressed value. The mass transport rate is lowered by their ratio, 
which is the incoming Mach number. 
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Since the contraction speed in the outer portion of the dense core is well under c s , the reduction 
in M may be substantial. The accretion luminosity is similarly reduced from traditional values, 
alleviating the discrepancy between theory and observation in this regard. We further show that 
the region of infall spreads out at a speed that is close to the subsonic injection velocity in the 
core's exterior. This relatively slow expansion helps to explain why global, supersonic infall has 
been difficult to observe. 

In Section 2 below, we describe the physical problem to be solved, present the relevant equations 
in nondimensional form, and then give our method of solution. Section 3 presents numerical results. 
Finally, Section 4 compares our findings with those of other studies, and indicates future directions 
for extending this work. 



2. Formulation of the Problem 
2.1. Physical Assumptions 

Dense cores are supported against self-gravity not o nly by thermal pressure, but also by the in- 



terstellar magnetic field, which penetrates their interior (ICrutcherll 19991 ). The lim ited data from po 



lariza tion mapping indicate that the field has a relatively smooth spatial variation (|Rao. Girart. Marrone 



20111 ). as would be expected from the absence of significant MHD wave motion. From a theoretical 
perspective, the objects are essentially magnetostatic structures. They must evolve quasi-statically, 
however, since at some point they undergo gravitational collapse. This evolution is facilitated both 
by ambipolar diffusion, i.e., the slippage of neutral gas across field lines, and by the addition of 
fresh material from the core's surroundings. 

There is no reason why the dense core's accrual of gas should cease once that object has 
produced a central star. Prior to this event, the core had an extended period of slow evolution. 
Thus, the internal accretion rate onto the protostar has had ample time to adjust so as to match 
the external rate at which mass is continually added to the dense core. In other words, M is 
determined by conditions far from the star itself. This was also the case in the traditional picture 
of an isolated cloud undergoing collapse, where M is set by the global free-fall time. A new picture, 
more consistent with observations, is that a dense core harboring a stellar source is a temporary 
repository for material that drifts in from the environment on its way to the protostar0 

To begin exploring the new picture qualitatively, we adopt a very simple model. We ignore 
magnetic forces entirely, and treat the dense core as a region of steady, isothermal flow surrounding 
the protostar. Recent observations by the Herschel satellite have shown vividly how dense cores 



The near equality of internal and external accretion rates does not hold if the external gas is injected supersoni- 
cally. In that case, simulations find that the inner rate underg oes a sharp, transient b urst just after star formation, 
before ultimately equilibrating to the outer rate (see Fig. 13 of iGong fc Ostriker 20091 ) . 
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are nested inside pc-scale filaments (jMen'shchikov et al.ll2010l ). Hence, the inward flow of gas is 
assuredly anisotropic. We ignore this complication as well, and assume spherical symmetry. Every 
fluid element experiences both the gravitational force from the protostar a nd tha t from interior 
cloud gas. Our problem is thus identical to the classic inflow calculation of iBondil (|1952l ). except 
that we include self- gravity of the gas. 

Like Bondi, we apply the inner boundary condition that this material is in free fall onto the 
star . However, we no longer demand that the density approach a fixed value at spatial infinity (as 
did IChialll978l ). Instead, we require that the velocity reach some imposed, subsonic value. Note 
again the critical difference from the many numerical collapse simulations, where the infall velocity 
was set to zero at some radius, effectively isolating the dense core from its surroundings. In effect, 
our imposition of a finite, bounding velocity means that we are treating protostellar accretion as a 
flow problem , rather than a collapse, which would necessitate a time-dependent treatment (see the 
discussion in IWandell Il984l ) . 



In reality, of course, the subsonic flow begins at a finite boundary. By taking this boundary 
to infinity, we are focusing, again for simplicity, on the earliest phase of protostellar accretion, 
when the mass of the central object is much less than that of the parent dense corej^l Since 
obser yations indicate that the final protostar mass is no m ore than 20 to 30 percent that of the 
core (jAlves. Lombardi. &: Ladal 120071 : iRathborne et al.1 [2009j) , this early epoch spans a substantial 
fraction of the protostar's total lifetime. 



2.2. Quasi-Steady Flow 

In spherical symmetry, momentum conservation in an isothermal, steady-state flow reads 



du G(M* +M r ) c 2 dp 
dr r 2 p dr 



u ~ = 2 7 — • ( 2 ) 



Here, u is the fluid velocity, taken to be positive for inward motion. The quantities M* and M r 
are, respectively, the mass of the central protostar and that of the cloud gas interior to radius r. 
Finally, p is the mass density of the gas and c s the isothermal sound speed. The cloud mass M r 
and density p are related through 

— - — = 4 7rr p . (3) 
dr 

Since the flow is steady, the same mass per unit time crosses every spherical shell. Denoting this 
mass transport rate by M, we have 

M = Airr 2 pu . (4) 



2 Specifically, we are requiring that M* <C rtCg/G, where rt, is the dense core boundary. The dense core itself is 
marginally unstable gravitationally, so the righthand side of this inequality is close to the core mass. 
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The inner boundary condition is that the velocity approach free fall onto the star: 

lim u = a/—— — * • (5) 

r- H) V r 

Equivalently, M r must approach zero in this limit. We set our outer boundary condition at spatial 
infinity and demand that 

lim u(r) = Uoo , (6) 
r— >oo 

where is some fixed, subsonic velocity that is consistent with the spectroscopic observations. 

Equation (4) incorporates our assumption that the inner accretion rate onto the star matches 
the outer rate imposed externally. In more detail, the match will not be exact, and the difference 
drives secular changes in the gas density and total core mass that only a full evolutionary calculation 
can track. Within the infall region, the assumption of steady state is justified as long as the crossing 
time through that region is brief compared to the time for M* to increase significantly. After 
completing the calculation, we will verify its self-consistency in this regard in Section 4 below. 



2.3. Nondimensionalization and Method of Solution 

It is convenient to recast our dynamical equations into nondimensional form. We adopt, as 
basic quantities, G, M*, and cs- The solution is then characterized by a dimensionless parameter, 
the Mach number associated with the inflow at the boundary: 

ft = u^/cs . (7) 

We define a nondimensional velocity as u = u/c s , and a nondimensional cloud mass as M r = M r /M* 
In addition, we define a fiducial radius, density, and mass transport rate, respectively, as 

GM * 

r Q = , (8) 

po = g^mI ' (9) 

Mo EE |, (10) 



and set r = r/r Q and p = p/p - 

Our dynamical equations become 



du 1 + M r 1 dp 

df f 2 p dr 



U ~ = i2 ' ( U ) 



= 4nr*p, (12) 
A = 47rr 2 pu, (13) 
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where A = M/M a . The boundary conditions are 

lim u = \l ^ , (14) 
f->0 V r 

lim u = p . (15) 

f — >oo 

From now on, we drop all tildes and alert the reader whenever we revert to dimensional variables. 
To solve this system, we first combine equations (12) and (13) into 

dM r A 



dr u 



(16) 



Next, we take the logarithmic derivative of equation (13) and use the result to eliminate the density 
gradient in equation (11): 

'„ _ i>) *. = _!+«: + h . (17) 



u J dr r 2 r 

The basic strategy is to solve equations (16) and (17) to obtain u(r) and M r (r). Equation (13) 
then provides an algebraic relation for p{r). 

Before describing the solution technique in more detail, we first derive a simple relation between 
A and (3. First note that equation (16) tells us that dM r /dr asymptotes to a finite value for large 
r. Thus, M r itself increases without limit. On the other hand, since u approaches a finite value, 
du/dr must vanish asymptotically. The entire lefthand side of equation (17) thus also goes to zero. 
The vanishing of the righthand side of this equation leads us to conclude that 

lim M r = 2r , (18) 

r— >oo 

confirming our previous claim that M r increases indefinitely. 
Equation (18) implies that 

lim ^ML = 2. (19) 

r— >oo dr 

We now apply this same, large-r limit to equation (16). After utilizing equation (15), we find the 
desired relation: 

A = 2 . (20) 

In dimensional language, this is 

M =(^)|. (21) 
Here we see explicitly the reduction factor modifying the conventional mass accretion rate. 

Supplied with this result, we proceed to solve the coupled equations (16) and (17) as follows. 
We start at an r- value inside the sonic point r s , i.e., the radius where u = 1. At this interior point 
rj, we guess both u and M r . Integrating outward, we find that u(r) diverges upward or downward as 
we approach r s , depending on the guessed n(rj). Using a bifurcation technique, we find that value 
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of u(rj) which extends the solution farthest out. We then integrate inward from rj, and similarly 
find that the velocity either climbs above or falls below the free-fall value given by equation (14). 
A similar bifurcation technique then yields the best-fit M r at r = r{. 

Having satisfied the inner boundary condition and established the location of the sonic point r s , 
our next task is to cross this point. We apply L'Hopital's rule to equation (17) and use equation (16) 
to find 

du Jl- A/2 

-T- = — , at r = r s . 22 

ar r 

This relation, together with the numerically established M r (r s ), allows us to integrate outward. As 
long as A obeys equation (20), the velocity u automatically approaches B at large r. 



3. Numerical Results 

Figure 1 displays the velocity profile u(r) obtained in the manner indicated. Here we have 
set the incoming Mach number equal to the representative value f3 = 0.2. Close to the origin, the 
velocity diverges as gas freely falls onto the protostar. Farther out, u first dips below B and then 
asymptotically approaches it through a series of decaying oscillations. 

The sonic point is found to be r s = 0.544, a slight increase over the Bondi value of 0.5. The 
dimensional sonic radius is then 

r s = \r s c s t (23a) 
= 2 r s /3 c s t . (23b) 



Thus, the point moves out with speed 



^ = 2f s f3c s . (24) 



If we identify the infall region as that volu me i n terior to r s , then this region expands at 0.22 c s . 



Note again the contrast with the model of IShul (119771 ). where the equivalent speed is exactly c s . 



Note also that the value of 2 f s is close to unity and insensitive to /3, rising only to 1.22 for (3 = 0.4. 
In summary, the infall region spreads out at a speed close to /3c s , at least during this early epoch. 

Figure 2 shows the density profile p{r). Inside r s , the profile is that associated with freely 
falling gas. From equations (13), (14), and (20), this limiting profile is 

r < r s . (25) 



r 2^/2vrr 3 / 2 

The dimensional equivalent is 



^ ~ (26) 



c s t 4vrGr 3 / 2 
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where we have also used M* = Mt. At early times, the density declines throughout the free-fall 
region, while the region itself expands. 

Outside of r s , equations (13), (15) and (20) tell us that 

P = r > r s . (27) 



which is dimensionally equivalent to 



c2 



' " • (28) 

This is the profile of the singular isothermal sphere, as we expect for a self-gravi tating gas cloud 



that contains a star of relatively small mass and is close to hydrostatic equilibrium (jStahler h Palla 



20041 . Chapter 9). 



4. Discussion 

Having solved the problem as posed, our next task is to check that the assumption of steady- 
state flow is reasonable. The outer portion of the cloud evolves slowly with time, as a structure 
that is never far from hydrostatic balance. It is the interior region that concerns us. Again, we 
need to make sure that the crossing time from r s to the origin is less than M*/M, which in turn is 
close to the evolutionary time tl§ Thus, we require 

• dr , . 

— < t . (29) 
u 

But u > c s in this region, so that the integral is less than r s /c s . Since the infall region expands 
subsonically, this ratio is indeed less than t. 

We stress the need for a fully time-dependent calculation to verify our main results. As- 
suming there is n o qua litative change, it is instructive to compare our findings with those of 



Gong &; Ostrikerl (120091 ). Using direct simulation, these authors built up a dense core from a 
spherically converging, supersonic flow. The accretion shock bounding the core first expands and 
then contracts supersonically before the core itself undergoes interior collapse. When the proto- 
star first appears, the velocity profile throughout the dense core is uniform and supersonic, with 
u ~ 3c s . A rarefaction wave erupts from the origin, but now expands supersonically. The wave 
quickly overtakes the accretion shock a nd effectively destroys i t, so that the converging inflow 



thereafter directly impacts the protostar. iGong fc Ostrikerl (|201ll ) repeated the simulation using a 



planar, external flow and found substantially the same results. 



3 The time M*/M would not be close to t if M varied rapidly in time, as it does in simulations of isolated, collapsing 
clouds. However, we do not expect such rapid variation when a subsonic, external flow is imposed from the start. 



- 10 - 



We observe once more the close relationship between the velocity of r s and that of the gas 
introduced externally. To see this connection in yet another light, note that r s is located roughly 
where the free-fall velocity onto the star matches the sound speed: 



2GM» 



(30) 



If we let M* = Mt and use equation (21), we find that 



Ur>n t. 



The fact that the outer regions of dense cores are contracting subsonically both before and 
after the appearance of a protostar appears to rule out a purely hydrodynamic origin for these ob- 
jects. The many researchers who have simulated cluster formation in turbulent molecular gas have 
effectively sidestepped this issue. In these studies, random velocity fluctuations are applied to a rep- 
resentative cube of molecular gas. Once these fluctations thorou ghly stir up the gas, self-gravity is 



switched on and overdense structures promptly collapse (see, e.g.j Ballestc ros-Paredes. H artm ann. k, Vazquez-Sema 



199fll : iHeitsch. Mac Low, fc KlesserJhoOll : badoan fc NordlundlboodJoffner. Klein, fc McKeelbooal ) 



More realistically, a dense core is self-gravitating and gains mass prior to this collapse, and a 
force in addition to the thermal pressure gradient must support it during this early epoch. This 
force is already known - the pressure gradient associated with the interstellar magnetic field. The 
field may not play a key role in t he collapse phase studied here, alt hough it apparently impacts the 
formation of circumstellar disks (ILi. Krasnoplolsky. Shangll201ll ) . In any case, the magnetic field 
is surely key in the early buildup of dense cores within their parent cloud filaments, and must also 
mediate the transition between the supersonic exterior of each core and its supersonic exterior. 



We are grateful to Paola Caselli, Neal Evans, and Phil Myers for sharing their expertise on 
dense cores and the subtle ways in which observations are revealing their dynamics. SWS was 
partially supported by NSF Grant AST-0908573. 
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Fig. 1. — Velocity profile in the quasi-steady flow, for the representative case (3 = 0.2. Both the 
velocity and radius are in the nondimensional units defined in the text. The velocity asymptotically 
approaches the imposed, outer value, shown as the dashed, horizontal line. 
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Fig. 2. — Density profile in the quasi-steady flow, for f3 = 0.2. The density and radius are both 
nondimensional. The dashed curve is a density profile that varies as r -2 . 



